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ABSTRACT 



^ . We review a recent attempt to deal with non-perturbative features of QCD by analytical 
O I means, using a manifestly gauge invariant, canonical approach. 

in 



INTRODUCTION 



in 

^ . Suppose you would like to solve the elementary quantum mechanical problem of a 

'- ^ ', particle in a central potential, but under the condition that only s-waves are "physical 

' ' states" 



L\phys} = . (2) 

This is a typical example of a constrained system, formulated in terms of redundant 
variables. Here, the constraint can easily be resolved by transforming to polar coordi- 
nates, 

f 1 \ 

\- — -^ + V{r)ju{r) = Eu{r), u{r) = r^{r) . (3) 

If we denote r again by x, the result is suggestive of the "axial gauge" y = z = (Note 
however that the condition x > and the boundary condition m(0) = for the radial 
wave function are remnants of the transition to curvilinear coordinates.) Rotational 
symmetry is guaranteed, irrespective of any further approximations to dynamics, since 
we are using the scalar variable r. 

In the case of gauge theories, we face a similar situation since we start out with re- 
dundant variables. Can we find the analogue of polar coordinates in QED or QCD, of 
course with respect to local gauge symmetry rather than rotational symmetry? Let us 
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consider QED first (Lenz et al, 1994a). The canonical formulation is most straightfor- 
ward in the Weyl gauge {Aq = 0) since Aq has no conjugate momentum. All 3 spatial 
components of A are then quantized. The Gauss law as constraint on the physical 
states accounts for the fact that only 2/3 of the variables (the transverse photons) are 
physical. The Hamiltonian in the Weyl gauge reads 



/" r 1 

H = J d^x ip^ (-iaD + (3m) i) + -{e"^ + 



(4) 



with 

D = y-ieA, E(x) = --^— , B = yxA. (5) 

« 5A{x) 

H is invariant under local, time independent gauge transformations 

A^A + WP, ^^e*'>, (6) 

generated by the Gauss law operator 

G{x) = -WE + ep , [G(f), H] = Q . (7) 

Physical states are defined through the constraint 

G{x)\phys) = Q . (8) 

Here, unlike in the above toy model, the choice of gauge invariant variables is far from 
unique, as could have been guessed from the proliferation of "gauge choices" in the 
literature. To see how one can resolve the constraint, write down the Gauss law in the 
Schrodinger representation, 

(vii + ep^vl/[l,V.] = 0. (9) 

This (functional) first-order differential equation can be solved by a "plane wave" 
ansatz, 

^[1, V^] = exp 1 d?xA^Vep^ <l>[^, V^] := f/^<f V'] • (10) 

Since U as defined in Eq. (|iy) is a unitary operator, the Schrodinger equation can be 
recast into the form 

= E'^ UHU^^ = E^ . (11) 

This is the equation corresponding to the radial Schrodinger equation (|]) in the above 
example, provided we can show that UHW contains only the physical variables A^j., ijj. 
Indeed, since the transformed Gauss operator 

UGix)U^ = V--^ (12) 

commutes with UHU\ it is clear that the unitarily transformed Hamiltonian can- 
not contain the longitudinal part of A any more. A simple calculation confirms this 
expectation, yielding 



UHU^ = J d^x ip^ [-iaDtr + /?m) + ^ (^El + - e^P^P 



(13) 
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We have thus rederived the standard Coulomb gauge Hamiltonian, including the famil- 
iar static Coulomb potential. The Gauss law constraint is resolved, and all redundant 
variables have been eliminated; in this quantum mechanical scheme of gauge fixing, U 
is denoted as "unitary gauge fixing transformation" (Lenz et ai, 1994a). In QED, the 
Coulomb gauge is clearly singled out on physics grounds, since it is only in this gauge 
that static sources decouple from the radiation field. However, many other gauges 
are conceivable; thus, for instance, the Gauss law can alternatively be resolved by the 
ansatz 

^[1, ^] = exp |-i J d'xA^^^ {-V±E^ + ep) | $[1^, ^j] . (14) 

This leads to the axial gauge {A3 — 0) which, however, is less convenient from the 
point of view of atomic or molecular physics, since even static charges "radiate" . 
In QCD, it is much less obvious whether such a procedure can be carried through and 
what the preferred gauge choice is (Lenz et ai, 1994b). The Weyl gauge Hamiltonian 
now reads 

j d^x [i)^ (-idD + /3m) V' + tr (e'^ + ^^)] , (15) 
with the chromo-electric and -magnetic fields (a = 1, N"^ — 1) 

E" = ---^r , = V X A" + -gr^^A^ X A" . (16) 

t dA"- 2 

The Gauss law is non-linear, 

(-DE + gp) \phys) = , {DEf = V^" + gP^^A^E" , (17) 

hence it is much more difficult to resolve than the abelian one. If we want to implement 
the Coulomb gauge for instance, we have to solve an equation of the type {E^ — V0) 

DV(j)\phys) = ...\phys) (18) 

for (f). This is impossible (non-perturbatively) since one cannot identify analytically 
the zero-modes of the 2nd order partial differential operator DV (Gribov, 1978). By 
contrast, in order to implement the axial gauge, one only needs to solve a first order 
ordinary differential equation, 

DsE^lphys) = ...\phys) . (19) 

Since one has control over the zero modes of D^, it becomes possible to derive a 
Hamiltonian formulated exclusively in terms of unconstrained, physical variables in 
the axial gauge (Lenz et ai, 1994b). The whole procedure is quite involved, and I 
will make no attempt to go into any technical detail, nor even to show you the final 
Hamiltonian in full glory. Instead, in the following two sections, I will concentrate 
on one particularly instructive part of the Hamiltonian, trying to exhibit some basic 
physics for which the axial gauge is advantageous (Lenz et ai, 1995). 



CONFINEMENT 

Let us come back to the axial gauge in QED for a moment and point out one compli- 
cation which we have ignored so far. We work in a finite box with periodic boundary 
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conditions (i.e., on a "torus"). In that case ^3 = is not a legitimate gauge, simply 
because the 2-dimensional field a3{x±) = Jq dzAs{x±, z) is gauge invariant. Physically, 
it describes transverse photons polarized in the 3-direction and propagating in the (1,2) 
plane. In QED, this can easily be cured: Use the gauge d^A^ = (rather than A3 = 0) 
by retaining a3{x±) and eliminate some other variables instead (Lenz et al, 1994a). 
In QCD, the corresponding gauge invariant quantities formed exclusively out of A3 are 
the eigenvalues of the spatial Wilson loop winding around the torus ( "Wilson line" ) , 

p^igJ^dzAs ^ Y^igaaLy^ ^20) 

(as: diagonal matrix). Here, the residual 2-dimensional variables cause much more 
trouble than in QED. In fact, most of the work needed to resolve the non-abelian 
Gauss law has to do with these {Nc — 1) lower dimensional, color neutral fields, and all 
the conspicuous non-perturbative features displayed by the gauge fixed Hamiltonian are 
somehow related to them (Lenz et al, 1994b). Among these features, most noteworthy 
is a Jacobian, reflecting the transition from Lie algebra (A3) to group elements (spatial 
Wilson lines) in the process of gauge fixing. 

Is this a purely technical matter, needed to properly define the theory in the infra-red, 
or is there some real physics associated with the 03? By construction, it is clear that the 
2-dimensional variables 03 are those whose dynamics has been maximally simplified by 
the choice of the axial gauge; as such, they may teach us something about the dynamics 
of a whole class of variables for which they are representative (but which are not as 
simply described in our "coordinate frame"). In order to exhibit their physics content, 
let us do a very drastic approximation - truncate the axial gauge Hamiltonian by 
keeping only the terms in 03 and the corresponding conjugate momenta (Lenz et ai, 
1995). Surprisingly, even this primitive version of the Hamiltonian still exhibits very 
interesting differences between QED and QCD, to which we now turn. 
Formally, the truncated Hamiltonian for both QED and SU(2) Yang Mills theory re- 
duces to the following 2-dimensional expression. 



/r 1 L ^ 

(fx — e^es + -(V^as) 



(21) 



Here, the electric field energy is found to be 



t 1 5 S 



J = sin' ( ) , (23) 



with the SU(2) Haar measure 



whereas J = 1 in the case of QED. This form of the kinetic energy requires some 
ultra-violet regularization, for which we choose a transverse lattice (lattice spacing ^, 
fundamental lattice vectors 6). (For QED, this would not be necessary but we do 
it anyway for ease of comparison.) Next, a "radial" wavefunctional is introduced as 
^ = "n/J^, and the kinetic energy is reduced to standard form like in the above quantum 
mechanical example, with a corresponding boundary condition (^ = whenever J = 
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0) . In terms of the rescaled variable (/? = the "radial" Hamiltonian then 
becomes 

h^he + hm^-J^Y.Q^ + ^T.i'Pr+S-Vr) ■ (24) 

Let us compare the abelian and non-abelian cases: 

1) QED: 

Since the Jacobian is trivial, the Hamiltonian is quadratic and can be diagonalised by 
a standard lattice Fourier transform (i^^ = Zlg e^''*^^/^^^, with N — L/l the number 
of lattice sites in each transverse direction.) As a result of the magnetic coupling, 
collective excitations appear, with the dispersion relation 

-i-|i:-^(f)-(^)^ (1*1 (25) 

In the limit L — > oo, they become just ordinary, massless photons. The ground state 
wave functional is Gaussian, 



k 



I 4£^ \ 
*-nexp(-^c.,-40,-) , (26) 



and the virial theorem ensures that fluctuations of electric and magnetic fields are 
equal, 

(^2 _ ^ Q ^ (^27) 

Hence the axial gauge is not a bad choice at all for pure QED: We have reduced 
the problem of solving a 3-dimensional vector field theory to that of a 2-dimensional 
scalar field theory, by an appropriate choice of coordinates. Although the Coulomb 
gauge is preferred for static charges, the axial gauge is well suited for studying the 
elementary gauge field excitations - a transverse field is naturally described in cartesian 
coordinates, the 3-axis pointing into the direction of the polarization vector. 

ii) QCD: 

Since ^ = whenever J = 0, we can restrict the variables </?f- to the interval [0, tt]. For 
large L, h^, dominates over /i^, and we are left with the simple quantum mechanics 
of infinite square well potentials at each site, totally decoupled from each other. The 
ground state is now uncorrelated in coordinate space, 

^0 ~ n Y - VP ■ (28) 

(Note that this corresponds to the "full" wavefunctional \&o =const., Eq = 0.) Since \l/o 
is an eigenstate of the electric field operator with vanishing eigenvalue, we find trivially 
an exact "dual Meissner effect" 

(E^) = . (29) 
The vacuum has a magnetic condensate reminiscent of the QCD vacuum, 

{E^ -B^) <0 . (30) 
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Most importantly, there are no such excitations as "plane wave gluons" ; the gap to the 
first excited state is 

Ai?=^-oo (L^oo), (31) 

so that these gluonic degrees of freedom get frozen in the thermodynamic limit. The 
magnetic contribution to the ground state energy can be evaluated perturbatively; one 
finds 

mm\0) = ^ - ij • (32) 

In contrast to the QED case, the vacuum has a precise value of the electric field (namely 
zero) and therefore the largest possible fluctuations in A ("stochastic" vacuum). The 
excitation gap ( pT]) can be understood by comparing it to lattice gauge theory: If 
one retains only ^3 and insists on Gauss's law, the only gauge invariant excitations 
are electric flux lines in the 3-direction around the torus. Their energy in the strong 
coupling limit of Hamiltonian lattice gauge theory (Kogut and Susskind, 1975) is 

(33) 

We thus recover the strong coupling string tension. Nevertheless, our approach is quite 
different from the standard lattice where one gets the same type of flux quantization 
and string tension also in QED, in the strong coupling limit (Rothe, 1992). Here, 
we obtain a qualitatively different behaviour in the two cases, the difference being 
exclusively due to the Jacobian. 

We have seen that the simple degrees of freedom of the axial gauge are useful in order 
to study the existence or non-existence of plane wave gauge bosons. We find strong 
indications that QCD does not admit chromoelectric fields with constant polarisation 
vector over large distances. So far, we cannot say anything about the scale involved - 
what is large? This would clearly require taking into account the other gluonic degrees 
of freedom as well and going through some renormalization procedure. Nevertheless, 
we can get some information about the relevant scale by indirect methods, using lattice 
results as input. This is important in order to further assess the possible usefulness of 
the axial gauge. 

DECONFINEMENT 

It is generally believed that QCD undergoes a deconfining phase transition to a quark 
gluon plasma at a temperature Tc ~ 150 — 200 MeV (Miiller, 1985). Can the axial 
gauge Hamiltonian formulation of QCD contribute anything useful to this issue? At 
first sight, the prospects look rather somber: Finite temperature field theory means 
compactification of the time direction, so that the Weyl gauge = is no longer 
allowed (for the same reason that the axial gauge A3 = is not allowed on the torus). 
Moreover, there is evidence from lattice gauge calculations that the spatial Wilson 
loops show area law behaviour even above (see e.g. Karsch, 1994). This seems 
to indicate that the "simple" variables 03 are not particularly sensitive to the phase 
transition, and that consequently we would have no advantage over other approaches 
by using the axial gauge. 
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Fortunately, nature provides us with a very elegant way out of these problems. As 

a matter of fact, one can study thermodynamic properties like the deconfining phase 
transition by working strictly at T = 0, but in a spatial box contracted in one direction 
(say, the 3-direction) to 

L3 = /? = 1/T . (34) 

This result is intimately connected to Lorentz invariance (or, rather, Euclidean 0(4) 
invariance) and would not hold in non-relativistic field theories. Since it is rather 
unfamiliar, let me explain it with a simple example: Consider the partition function 

J D[0]e-/"''^^^(^) (35) 

for a generic field theory in a finite box at finite temperature in 2 cases: 

/) Li,L2,L3>/3, 

//) /?,Li,L2>L3 . (36) 




(I) 



(11) 



Case (I) corresponds to a hot system in a large, isotropic box, case (II) to a cold system 
in a box contracted along the 3-direction. Covariance of the Euchdean action relates 
these 2 situations which differ by the interchange of 2 coordinates, xs and the Euclidean 
time T. Using the standard relations 



one finds for instance 



dp 



P — — ^ — in Z , 



II 



(37) 



(38) 



and vice versa. These symmetry relations have surprising and powerful consequences. 
By way of example, consider the Casimir effect for a massless scalar field with periodic 
boundary conditions (Toms, 1980). If the system is enclosed between two plates at 
distance L3 = d, there is an attractive force corresponding to a negative pressure 



'30d'i ■ 



(39) 



7 



On the other hand, the energy density for an ultra-relativistic ideal gas of scalar par- 
ticles at finite temperature is given by the Stefan Boltzmann law 



^ = ^T^ (40) 

These two laws of seemingly unrelated parts of physics are indeed mapped onto each 
other by the substitution d <-> 1/T. This observation - which has been exploited 
several times in the literature (cf. e.g. Toms, 1980; Koch et al., 1992) - opens a new 
perspective for studying thermodynamic properties of field theories in a technically and 
conceptually simpler way. All one has to do is to study the ground state in a different 
geometry. 

Once this is understood, it is not difficult to identify other relations between at first 
sight unrelated physical observables. One particularly amusing example is the fol- 
lowing: In (Lenz et al, 1991), the Schwinger model was studied as a function of a 
parameter which measures how far off the light-cone the quantization surface is. The 
fermion condensate was evaluated analytically as a function of this parameter. It was 
also pointed out that the same function can be interpreted as dependence of the con- 
densate on the size of the box. Independently, the finite temperature Schwinger model 
was investigated in (Sachs and Wipf, 1992), and again the fermion condensate evalu- 
ated in closed form as function of T. The two formulae (Eq. (3.105) of (Lenz et al, 
1991) and Eq. (5.10) of (Sachs and Wipf, 1992)) agree exactly if one identifies the 
corresponding variables, 

< = lif3m,r , (41) 

a fact which seems to have been overlooked so far. Similarly, it is tantalizing to rein- 
terpret corresponding findings for large Nc QCD2 in (Lenz et al, 1991) as evidence for 
a chiral symmetry restoring phase transition at finite temperature, in contradiction to 

the common lore (McLerran and Sen, 1985; Ming Li, 1986). 

Let us now return to axial gauge QED and QCD. We start with QED and go one step 
beyond the truncation of H, taking into account perturbatively the coupling of to 
the charged fermions. The simplest way to do this is to compute the effective potential 
for the zero mode of 03 by evaluating the energy density of the Dirac sea in a constant 
background potential 03. In view of the application to finite temperature field theory, 
we have to require anti-periodic boundary conditions for the fermions, 

V'(f±,L) = -V'(f±,0) . (42) 

We can gauge away a constant as, provided we change these boundary conditions into 
quasi-periodic ones, 

V'(f±,L) = -e-"^^V(^±,0) . (43) 
This yields the following discretization for the 3-component of the fermion momenta, 

P3 = ^(2n + l) -eag . (44) 

The effective potential is then given by 

Uesias) = -2C/(c) (45) 
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with U{c) the (heat-kernel regularized) sum over single particle energies, 

U{c) = / ^ E E{p^^ ^ - c)e-^^(^^— ) , (46) 



The variable c is defined as 



(2vr)^ 

eLa^ 1 
~ 27r ~ 2 



(47) 



and 



E(p-1, u) = ^pI + (27rz//L)2 . (4J 



Performing the integral in (^) and using the generating function for Bernoulli poly- 
nomials, one finds 

1 q2 1 oo 

^ ^ 27rL dX' X ^^too 

27r2 

54(c) (49) 



3L4 

(valid for |c| < 1, and to be continued periodically outside this interval). We have 
dropped a c- independent, infinite constant. The result for the effective potential is 
therefore 

47r2 

^^efr(a3) = -^^4(c) 

iiig, LU t; — 2 '-'J- "3 — 

Expanding t/gff to 2nd order around the minimum, we find 



f/cff has extrema at c = 0, |, 1, the minimum corresponding to c = | or 03 — ^j^. 



From this expression, replacing L by /5 = 1/T, we can read off the familiar free energy 
density of an ideal gas of massless fermions, 

as well as the "electric mass" of the photon (Kapusta, 1989), 

^ei = \e^T^ , (53) 

in spite of the fact that we have only been dealing with ground state properties. This 
gives a first impression of the potential use of the axial gauge, if at the same time one 
reinterprets the contracted box in terms of finite temperature field theory. 
Let us now turn to the corresponding computations for QCD. Here, if one interchanges 
3- and time directions, the simple variables in the axial gauge become the eigenval- 
ues of the thermal Wilson lines, the standard order parameters for the confinement- 
deconfinement transition. This is obviously a welcome feature which means that we 
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have chosen a useful set of variables and are in a good position to derive an effective 
theory for this order parameter in the spirit of Landau Ginzburg theory. 
Ignoring quarks (which could be taken into account easily as well in the present ap- 
proximation), we evaluate the zero-point energy of the "perpendicular" gluons A\^Ai 
in the presence of a constant background field 03. We treat the as free fields, except 
for the minimal coupling to 03. Again, if 03 is constant, its effect is equivalent to chang- 
ing the boundary conditions in 3-direction from periodic to quasi-periodic ones. The 
expression for the vacuum energy density can be written down in complete analogy to 
the QED example. We consider SU(2) Yang-Mills theory and introduce the rescaled 
field variable 

c = |^a3 G [0, 1] . (54) 

The vacuum energy density of A^ will serve as effective potential for 03. It can be 
decomposed as 

f/eff = 2f/(c) + f/(0) , (55) 

where the two different terms come from the (color) charged and neutral components 
of A^, respectively. The function t/(c) is exactly the same as the one given in Eq. (^6|). 
This yields for ?7efr, Eq. ([55|), the form 

t^efr = ^fc^(l-c)^-T^) (56) 



3L4 V 20 

(here again, this function has to be continued periodically outside of the interval < 
c < 1). The value of t/gff in the minima reproduces the free energy density of an ideal 
gluon gas, if we replace L by = 1/T, 

f/eff|c=0 = -yI^ - -^T^(iVe' - 1) (iVe = 2). (57) 

Using the same substitution, we reproduce the "one-loop effective potential" for the [x- 
independent) thermal Wilson line, which has been discussed extensively in the literature 
(first for SU(2) by Weiss, 1981). A look at the original reference shows that our 
derivation is significantly simpler. 

In the standard approach to finite temperature QCD, one infers the gluon electric mass 
from the 2nd derivative of in a minimum. 



m, 



, 1 0^2 



„ 3L2 3 

c=0 



N,g'T' (58) 



From this point of view, it would seem that QED and QCD behave quite similarly 
indeed. However, so far, we have only discussed the effective potentials and ignored 
the difference in the kinetic energy of as, which was crucial to "confine" plane wave 
gluons. In our approach, we see no justification for disregarding these effects, which 
are overwhelming at low temperature. 

Let us try to understand at least qualitatively the effect of the Jacobian. Following 
(Polchinski, 1992), we first assume that the effective potential derived for a constant 
as can be taken over for x^-dependent 03 as well, with the same functional form. The 
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effective potential to be added to hg + hm, Eq- (k4|), is tlien (up to an irrelevant constant 
term) 

Ucs = vl (tt - Vff ■ (59) 

r 

Secondly, we model the effect of the Jacobian by approximating the infinite square 
well potential by a harmonic oscillator potential, chosen such as to reproduce the gap 
between ground state and first excited state, eq. (0), 

«conf = ifr - TT/2f . (60) 

The resulting potential Utot = u^s + Wconf now exhibits the correct Z2 "center symme- 
try" , whereas without the confining potential it has a discrete translational symmetry 
both in QED and QCD. With a quartic potential and a quadratic one (possibly chang- 
ing sign). 



Mtot = ( - 1^ j Y i'^r - vr/2)^ + Y iv^r " vr/2)^ + const. , (61) 



we are now in the standard situation for a second order phase transition and can 
determine the gluon mass below and above the critical temperature. It is convenient 



to introduce 2 masses, the electric gluon mass of Eq. (|58|) and the "confining mass" of 
Eq. (^, 

-^el=^, -conf=(^^j ■ (62) 

The critical temperature (using L = 1/T) where the quadratic term in Utot vanishes is 
then determined by the condition m^^^f — m1^/2 = 0, or 

The expressions for the gluon mass in the confined and deconfined phase are = 
"^conf ~ and m^g = m^j — 2m^Q^£, respectively. These results are by no means 

realistic, but are reminiscent of discussions of the phase transition in the strong cou- 
pling limit (Gross, 1983). Nevertheless, it is interesting that the crudest conceivable 
approximation to the gauge fixed Hamiltonian has already the potential of describing 
both confinement and deconfinement of gluons. Much more work is needed in order to 
properly understand how the confining effects (related to the Jacobian) are overcome 
at high temperature and how the well established perturbative features of QCD can 
be recovered in the present framework. 

We finish with the following simple observation, which is again related to the question of 
length scales. If we translate the critical temperature as determined from lattice gauge 
calculations (Karsch, 1994) into a critical length, we learn that the phase transition to 
the quark gluon plasma takes place if we contract space in 3-direction to 

Lc = 1/Tc ~ 1.0 - 1.3 fm . (64) 

Thus, the variables 03 should be qualitatively representative for a large class of waves 
indeed - all those whose polarization vector does not change over distances appreciably 
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larger than 1 fm. On the other hand, Lc also sets the scale for the length at which 
quark and gluon degrees of freedom become essential - a rather large value from the 
point of view of nuclear physics. 
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